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Comparison of structure functions measured in neutrino and charged lepton deep inelastic scattering can be used to test basic 
properties of parton distribution functions such as the validity of charge symmetry. Recent experiments indicate a substantial 
discrepancy between FX and F% in the region of small Bjorken-x. We discuss nuclear corrections and strange and anti strange quark 
effects and show that none of them can account for the observed discrepancy. These results suggest surprisingly large CSV effects 
in nucleon sea distributions. 



1 Introduction 

Comparison of structure functions measured in neutrino 
and charged lepton deep inelastic scattering can be used 
to test basic symmetry properties of parton distribution 
functions such as the validity of charge symmetrylil. 

Such comparisons are based on the interpretation of 
these structure functions in terms of parton distribution 
functions. In the quark-parton model the structure func- 
tions measured in neutrino, anti neutrino and chacged 
lepton DIS on an iso-scalar target, Nq, are given bya 

F% N ° (x) = x[Q(x) + 2s(x) - Su(x) - 5d(x)} 



2 °(x) = x[Q(x) + 2s(x) - 8d{x) - Su(x)} 

ZN 



r^\r) = ^x[Q{x) + \{s{x)+s{x)) 



4(5d(x) + 8d{x)) (Su(x) + Su(x)) . 



(1) 



Here, we introduce the notation Q(x) = u(x) + u{x) + 
d(x) + d(x) and express everything in terms of parton 
distributions functions in the proton and charge sym- 
metry violating distributions which are defined as the 
differences between the up (down) quark distribution in 
the proton and the down (up) quark distribution in the 
neutron 



6u(x) 
5d(x) 



u p (x)-d n (x) 
d p (x)-u n {x). 



(2) 



A deviation R c (x) ^ 1, at any value of x, must arise 
either from CSV effects or from s(x) ^ s(x). 

Recent experimental measurements allow a precise 
comparison between F»(x,Q 2 ) and FJf(x, Q 2 ). In Figjl] 
we show the 'Jcharge ratio" calculated by using the CCFR 
and NMC □ structure functions for F 2 V and F$ , re- 
spectively. The structure functions are integrated over 
Q 2 = 3.2 GeV 2 in the region of overlap of the two ex- 
periments. Since the CCFR structure function has been 
measured on an iron target nuclear corrections (nuclear 
EMC effect, shadowing and anti shadowing) have to be 
applied to the data. In Fig. 1, the open triangles repre- 
sent the result without nuclear corrections. The open cir- 
cles are calculated by using nuclear corrections obtained 
from charged lepton DIS. While, in the region of interme- 
diate values of Bjorken x (0.1 < x < 0.4), the two struc- 
ture functions are in very good agreement, in the small 
x-region (x < 0.1), they differ by as much as 10-15% 
when nuclear corrections are taken into account. This 
discrepancy could be interpreted as evidence for charge 
symmetry violation. However, several corrections have 
to be applied to the data before any conclusions may be 
drawn. We see especially that the result is very sensitive 
to nuclear corrections. The CCFR Collaboration made 
a careful study of overall normalization, charm threshold 
and iso-scalar correction effects 13. Here, we re-examine 
nuclear corrections for neutrinos and discuss s(x) ^ s(x) 
effects before we turn to possible charge symmetry vio- 
lation. 



If charge symmetry is valid these terms are zero. 

A very useful quantity in the comparison of the struc- 
ture functions is the "charge ratio" 



Rc(x) 



F% No {x) 



18 F 2 



x[s(x) + s(x)}/6 



(3) 



2 Nuclear corrections 

Nuclear corrections for neutrinos are generally calculated 
using correction factors from charged lepton reactions at 
the same kinematic values. A priori, there is no reason 
that neutrino and charged lepton heavy target correc- 
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Figure 1: The "charge ratiofi R c vs. 2 calculated using CCFRQ 
data for neutrino and NMCu data for muon structure functions. 
Open triangles: no heavy target corrections; open circles: v data 
corrected for heavy target effects using corrections from charged 
lepton scattering; solid circles: v shadowing corrections calculated 
in the "two phase" model. Both statistical and systematic errors 
are shown. 



neutrino case, we expect that shadowing due to VMD 
in neutrino reactions is roughly half of that in charged 
lepton scattering. For large Q 2 -values, shadowing due to 
Pomeron exchange (which is of leading twist) becomes 
dominant, leading to identical (relative) shadowing in 
neutrino and charged lepton DIS. n 

We calculated shadowing corrections^ to the CCFR 
v structure function in the framework of this "two phase" 
model and used these corrections in calculating the 
charge ratio R c of Eq. |[ There are also nuclear effects 
in the Deuteron. However, because of the low density of 
the Deuteron, these are (relatively speaking) very small 
and have a negligible effect on the charge ratio. We inte- 
grated the structure functions above Q 2 = 3.2 GeV 2 in 
the overlapping kinematic region of the two experiments 
and used a parametrization of the nuclear corrections 
in charged lepton DIS to correct the data in the non- 
shadowing region. The result is shown as solid circles in 
Fig.0. At small x, careful consideration of neutrino shad- 
owing corrections decreases, but does not resolve, the 
low-x discrepancy between the CCFR and NMC data. 



tions should be identical, especially if such corrections 
depend strongly on the properties of the exchanged ob- 
ject (photon, W) used to probe the structure of the tar- 
get. Since this is the case for nuclear shadowing cor- 
rections in the small xb region for small to moderately 
large Q 2 -values we re-examined shadowing corrections to 
neutrino DIS focusing on the differences between neu- 
trino and charge lepton scattering and on effects due to 
the Q 2 -dependence of shadowing. We used a two phase 
model which has been successfully applied tp the descrip- 
tion of shadowing in charged lepton DIS u. In general- 
izing this approach to weak currents, subtle differences 
between shadowing in neutrino and charged lepton DIS 
arise because of the partial conservation of axial currents 
(PCAC) and the coupling of the weak current to both 
vector and axial vector mesons. PCAC requires that that 
the divergence of the axial current does not vanish but is 
proportional to the pion field for Q 2 = 0. This is Adler's 
theoremtJ, which relates the neutrino cross section to the 
pion cross section on the same target for Q 2 = 0. Thus, 
for low Q 2 w to 2 shadowing in neutrino scattering is 
determined by the absorption of pions on the target. 

For larger Q 2 -values the contributions of vector and 
axial vector mesons become important. The coupling of 
the weak current to the vector and axial vector mesons 
and that of the electro-magnetic current to vector mesons 
are related to each other by the 'Weinberg sum rule" 
f 2 + = f 2 x = 2/ 2 . Since the coupling of the vector (ax- 
ial vector) mesons to the weak current is twice as large 
as the coupling to the electro-magnetic current but the 
structure function is larger by a factor of ~ 18/5 in the 



3 Strange and anti strange quark effects 

Since the CCFR-Collaboration uses both neutrino and 
anti neutrino events in the structure function analysis the 
extracted structure function F§ CFR is a flux weighted 
average between v and v structure functions S 

F 2 CCFR (x, Q 2 ) = aFZ(x, Q 2 ) + (1 - a)F*(x, Q 2 ). (4) 

Here, we define the relative neutrino flux as a — 
^v/{^u + where <£>„ and <&p are the v and v fluxes, 
respectively. F§ CFR is equal to \{F%(x, Q 2 )+F%{x, Q 2 )] 
if a = i or if the two structure functions are equal, i.e. 
charge symmetry is valid and s(x) = s(x). The value of 
a depends on the energy of the incident neutrinos and 
anti neutrinos in the CCFR-experiment. In the relevant 
kinematic region, at small x which corresponds to large 
incident energies, it is « 0.83 so to a good approximation 
F 2 CCFR {x,Q 2 ) can be regarded as a neutrino structure 
function. 

The remaining small- a; discrepancy in the charge ra- 
tio is either from different strange quark distributions 
s(x) ^ s(x) or from charge symmetry violation. First, 
we examine the role played by the strange quark dis- 
tributions. Assuming charge symmetry, s(x) and s(x) 
are given by a linear combination of neutrino and muon 
structure functions, 

\f^ cfr {x, Q 2 ) -3F» mc (x, Q 2 ) = \ x [s(x) + s(x)\ 
+~(2a-l)x[s(x)-s(x)}. (5) 



2 




Figure 2: xs(x) (open circles) and xs(x) (solid circles) extracted 
by combining CCFR and NMC structure functions withs(x)' JA ' 
extracted from dimuon experiments, as given in Eqs. a and rj. Solid 
triangles: ^Ff CFR — 3F^ fMc . Stars: ozs(x) from a LO-analysis 
Q. Solid line: xs(x) from a NLO-analysis|j; dashed band indicates 
±lcr uncertainty. 



Under the assumption s(x) = s(x), this relation could be 
used to extract the strange quark distribution. However, 
as is well known, s(x) obtained in this way is inconsistent 
with results extracted from independent experiments. 

Opposite sign dimuon production in deep inelastic v 
and v scattering provides a direct determination of both 
six) and s(x). The CCFR Collaboration extracted s(x) 
and s(x) from aJeading order (LO)E_ next to leading order 
(NLO) analysis., of their dimuon data. While the strange 
and anti strange distributions were different in the LO 
analysis they were equal within experimental errors in 
the NLO analysis. 

To test the hypothesis that the low- a; discrepancy 
in the charge ratio of Eq. || could be accounted for by 
allowing s(x) ^ s(x) we combined the dimuon production 
data, averaged over v and v events, with the structure 
functions from neutrino and charged lepton scattering 
(Eq. |). Defining a' = N V /{N U + N p ), where N v = 
5,030, N„ = 1,060 (a' w 0.83) are respectively the J/ 
and v events from the dimuon production experiment □, 
the flux- weighted experimental distribution xsix) 11 ^ from 
dimuon production is 



xs^{x) = -i [s(x) + s(x)] 



+ -(2a' - l)x[s(x) - s(x)]. 



(6) 



Eqs. (|^) and @ form a pair of linear equations which 
can be solved for s(x) and s(x). We can simultaneously 
test the compatibility of the various experiments. 

Compatibility of the two experiments requires that 
physically acceptable solutions for ^x[s(x) + s(x)] and 



^x[s(x) — s(x)], satisfying both Eq. |5| and Eq. ^|, can be 
found. Clearly, solutions do not exist if |(2a— 1) = (2a' — 
1). (Note that the left-hand sides of Eq. || and Eq. H are 
different.) This gives the critical values a — 4 (3a' + 1) 
or a w 0.7 for a' « 0.83 and a = 0.8 for a' = 1. 

In Fig. |^ we show the results obtained for xs(x) (open 
circles) and xs(x) (solid circles) by solving the resulting 
linear equations, Eqs. || and ^, with a ~ a' « 0.83. Both 
the structure functions and dimuon data have been inte- 
grated over Q 2 > 3.2 GeV 2 in the overlapping kinemat- 
ical regions. In averaging the dimuo.n_.data we used the 
CTEQ4L parametrization for s fJ - fJ -(x)M. We see that the 
results are completely unphysical, since the equations re- 
quire s(x) < 0. Since a « 0.83 is closer to the critical 
value for a' — 1 than for a' — 0.83 the results is even 
more unphysical if we use only the neutrino dimuon data. 
Similarly, if we decrease a s(x) becomes more negative. 

In conclusion, our analysis strongly suggests that re- 
quiring charge symmetry, but allowing s(x) s(x), can- 
not resolve the discrepancy between F% CFR (x, Q 2 ) and 
F^ MC (x,Q 2 ). The experimental results are incompati- 
ble, even if s(x) is completely unconstrained. 



4 Charge symmetry violation 

We have seen that neither neutrino shadowing correc- 
tions nor allowing s(x) ^ s(x) removes the low- a: dis- 
crepancy. There remain the two possible explanations. 
Either one of the experimental structure functions (or 
s(x)) is incorrect at low x, or parton charge symmetry is 
violated in this region. Assuming the possibility of par- 
ton CSV, we can combine the dimuon data for s(x), (Eq. 
||), with Eq. H to obtain 



£ (x, Q z ) - 3i^2 (x, Q ) - xs^{x) 



\[\{2a-l)-{2a' -l)}x[s{x)-s{x) 
+— x [Sd(x) — Su(x)]. 



(7) 



This equation is valid at small x, where sea quark dis- 
tributions are much larger than valence quarks, so we 
make the simplest assumption, namely that 5q v {x) — 
Sq(x) — Sq(x) w 0. 

The left hand side of Eq. (?] and the coefficient in 
front of s-s are positive. Consequently, the smallest CSV 
effects will be obtained when s(x) = 0. Note that this 
violates the requirement that the net strangeness of the 
nucleons be zero. Thus, setting s(x) = gives the ab- 
solute lower limit on charge symmetry violation. In Fig. 
U we show the CSV effects needed to satisfy the exper- 
imental values in Eq. ^. The open circles are obtained 
when we set s(x) — 0, and the solid circles result from 
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Figure 3: Charge symmetry violating distributions extracted from 
the CCFR and NMC structure function data and the CCFR 
dimuon production data under the assumption that s(x) = s(x) 
(solid circles) and s(x) « (open circles) for a' = 0.83 and 
s(x) = s(x) (solid circles) and s(x) PS (open triangles) for a' = 1. 
(For the latter only statistical errors are shown.) s(x) at Q 2 rt 4 
GeV 2 obtained by the CCFR Collaboration in a NLO analysisEI is 
shown for comparison. 



setting s(x) — s(x). We also show the results we ob- 
tain by setting a' = 1 which corresponds to using only a 
subsample of the di-muon data containing only neutrino 
events. For s(x) = s(x), the result is the same. However, 
setting s(x) = gives us a higher lower limit on CSV 
than that for a'. The result is shown as open triangles 
in Fig. |. 

The magnitude of the extracted CSV is also sensi- 
tive to the strange quark distribution used in the anal- 
ysis. In Fig.[| we show the CSV obtained by using the 
LO (open circles) and NLO (solid triangles) CCFR dis- 
tributions and the CTEQ4L (solid circles) and CTEQ4D 
(solid rectangles) parametrizations. The uncertainty due 
to different parametrizations has been partly taken into 
account since the calculated errors already include the 
uncertainty of the dimuon measurement and most of the 
parametrizations lie within the experimental errors of the 
dimuon data (except of LO-CCFR s(x)). 

In conclusion, the CSV effect required to account 
for the low- a; NMC-CCFR discrepancy is extraordinarily 
large. It is roughly the same size as the strange quark 
distribution at small x. This CSV term is roughly 25% 
of the light sea quark distributions for x < 0.1, and the 
sign gives d p (x) > u n (x) and d n (x) > u p (x). 

5 Influence of CSV on other observables 

Clearly, CSV effects of this magnitude need further ex- 
perimental verification. It is hard to imagine how such 
large CSV effects are compatible with the high precision 



Figure 4: Uncertainty due to the parametrization used for the 
dimuon data on charge symmetry violation. Open circles: LO 
CCFR distribution, solid circles: CTEQ4L, solid rectangles: 
CTEQ4D, solid triangles: NLO CCFR distribution. Here, except 
for the most "critical" point, only statistical errors are shown. 



of charge symmetry measured at low energies. The level 
of CSV required is surprising, as it is at least two orde 
of magnitude larger than theoretical CSV estimates tilt 
Theoretical considerations suggest that 6d(x) « — Su(x) 
E£l; with this sign CSV effects also require large flavor 
symmetry violation. Since most of the observables are 
proportional to the sum of 5d(x) and 8u(x) rather than 
to their difference this large CSV could remain unob- 
served in many experiments. Here, we briefly discuss the 
effects of the extracted CSV on FSV. 

The Drell-Yan experiment measures the ratio of the 
dimuon cross sections of proton-Deuteron and proton- 
proton scattering. Since CSV is significant in the small 
x region, it is a good approximation to keep only the 
contributions to the Drell-Yan cross sections which come 
from the annihilation of quarks of the projectile and anti 
quarks of the target. In this approximation, the ratio 
R = cr pD /(2a pp ) is given by 



a pD 
2a pp 



1 h I d2 Sd2 ] -L Rl [1 I d 2 Su 2 1 
1 ^ ~r U 2 U2 ' 4 ^ ' fi 2 «2 1 



1 



ill §2 

4 u 2 



(8) 



Here, we introduced the notation Rx = and gi_2 = 
9(^1,2) for the quark distributions. Neglecting R\ for 
large xp, which corresponds to large x\, we have 



R = 



2<j pp 



1 



{1 



(rig - Sd 2 

U2 



(9) 



If charge symmetry is violated, the extracted quantitylll 
is not dijui but r-i = (dq — 5d2)/u2- Since Sd is positive 
the FSV ratio, ^2/^2, should be enhanced for small x 
relative to the measured value r%. The enhancement is 
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Figure 5: The ratio d/u extracted from the Drell-Yan data assum- 
ing the validity of charge symmetry (solid circles). If CS is violated 
this ratio corresponds to (d — Sd)/u. The result obtained by cor- 
recting for CSV is shown as open circles. The same for the ratio 
extracted from the difference of the proton and Deuteron structure 
functions is shown as solid and open triangles, with and without 
CSV, respectively. 



In conclusion, we have examined in detail the discrep- 
ancy at small x between the CCFR neutrino and NMC 
muon structure functions. The only way we can make 
these data compatible is by assuming charge symmetry 
violation in the sea quark distributions. The CSV am- 
plitudes necessary to obtain agreement with experiment 
are extremely large - at least two orders of magnitude 
greater than theoretical predictions of charge symmetry 
violation. If CSV effects of this magnitude are really 
present, then one must include charge symmetry violat- 
ing quark distributions in phenomenological models from 
the outset, and re-analyze the extraction of parton dis- 
tributions. 
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in the order of 25% in the small x region where CSV 
could be important. This is shown in Fig. 

The difference, d — u, can also be extracted from 
the difference between the proton and neutron structure 
functions measured by the the NMC Collaboration^ In 
this case we have 

-(u v -d v )-— (ff-J?) = {d-u)--(Sd+Sd)—(Su+6u). 
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(10) 

We can make the approximations Sq a Sq and Sd a — Su, 
and obtain 



2 K " d v ) - — (ff - F 2 ") a [(d- 5d) - «]. (11) 
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